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Abstract. From the observing campaigns of a number of 
hclioseismic telescope networks such as the Global Oscil- 
lation Network Group (GONG) and also from the Solar 
Heliospheric Observatory satellite (SoHO), helioseismol- 
ogists now have data on in excess of 20000 oscillation 
frequencies v of the modes known collectively as the 5 
minute oscillations. Depending on the order and degree 
of the modes these data are available either as individual 
frequencies of modes with known radial order n, spher- 
ical harmonic degree I, and order m, or as so-called 'a- 
cocmcients' which are obtained by fitting polynomials to 
ridges of power in the (/, v) diagram. Although there is a 
one-to-one relationship between these two types of data, 
the conversion between them is not straightforward which 
complicates comparison of helioseismic inference between 
the two types of data. Their relationship with each other 
and with the solar rotation rate is clarified in this paper. 
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1. Introduction 

For observations that are well resolved in space and in 
time the oscillations of the Sun can be projected uniquely 
onto its pulsation eigenmodes, which are products of func- 
tions of radius and of spherical harmonic functions. Each 
mode, and therefore each measured oscillation frequency, 
is uniquely identified by three numbers : the radial or- 
der n, and the degree / and the azimuthal order m of the 
spherical harmonic. For a non-rotating star the frequency 
is degenerate with respect to m. Rotation removes this 
degeneracy and the frequencies are split multiplets. Indi- 
vidual mode splittings D n i m (m > 0) are related to the 
mode frequencies by : 



2m 
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and to the solar rotation rate f2(r, 9) by 
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where r is the fractional radius and 9 the colatitude. The 
K n [ m are the mode kernels for rotation given by : 

K nlm (r,8) ee F?{r)G l r{6) + F^(r)G l 2 m (9) (3) 

where 

Ff ee p( r y [e n i{r)-2L- 1 U{r)rim{r)+ri l {r)) /I nl {A) 
and 



p(r)r 2 [L- 2 r, 2 m(r)] /Im 



(5) 



Here p is the density as a function of radius, L 2 = 1(1 + 1) 
and £ n i, 7] n i are the radial and horizontal components of 
the displacement eigenfunction in the non-rotating star, 
and 
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Further (using u = cos 9) 
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Here the P/™(m) are associated Legendre polynomials. 
This separation of K nlm into F? l ,F£ l and G[ m ,G l 2 m is 
not unique. Many other choices are possible but this choice 
is particularly useful for the purposes of this paper. The 
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same choice is made in Pijpers & Thompson (1996). One 
advantage of this choice is that : 



1 / 2^d 2 G ^ 1 " 



G l 2 m = 



Note that with the definitions (7), (8) of G l {^ the normal- 
ization is such that : 



J AuG\ m {u) = 1 = -J duG l 2 m (u) 

-l -l 
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or : 



(9) m=-i 
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and using the orthogonality of the V's then leads to : 

l i 2 

E vmmV<iXm) = a nls E (^ (0 (m)) (16) 



2. The a-coefficient fitting 

The a-coefficients are constructed by fitting polynomials 
in m to ridges of power in the (I, v) diagram. Details of 
their use can be found in eg. Schou et al. (1994, hereafter 
SCDT). The relationship between a-coefficients and mode 
frequencies is : 



Vnlm = E a nlsVf ] {m) 



s=0 



Here the polynomials V^p have degree s and are often nor- 
malized so that Vs l \l) = I- As noted in SCDT, who follow 
the projection proposed by Ritzwoller & Lavely (1991), 
these properties together with the orthogonality condi- 
tion : 



E Vi l \m)r { J } (m) = for s ^ s' 



(12) 



specify the polynomials completely. As noted in Ap- 
pendix A of SCDT the polynomials are usually con- 
structed in an iterative process. 

In practice it can be impossible, eg. for high values of 
I, to identify the individual peaks of power of every mode, 
in which case the fitting must proceed using a slightly 
different scheme than that described below. However at 
least formally the fitting always can be described as a 
least squares fit of polynomials to individual frequencies. 
In such a least-squares fit for the a-cocfficicnts the \ 2 ~ 
mcasure to be minimized is : 



X 2 = E \ v ni m -^a n i s V { i ] {m) 



(13) 
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Differentiation of Eq. (13) with respect to a n i s yields : 

0= E klm-Eanl^Hl^H (14) 



from which is obtained : 

a nls — , , 9 — 2_^i C lm Vnlr 



(17) 



where the are defined as 
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In the case of inversion for the determination of the solar 
rotation rate only the a-coefficients for odd s are of inter- 
est since these are a measure of the flow that is symmet- 
ric between northern and southern hemispheres. Since the 
coefficients satisfy cp_ s + 1} = -cf^ +1) (i.e. v£ +1 {-m) = 

~ ^2s+i( m ))> ^ is clear that linear coefficients 7^ +1 exist 
such that : 



a n l2s+l — E "fTs+lDnlm 



(19) 
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and they can be obtained from the polynomials by 
using : 



l2s+l 



m H l s+i( m ) 



(20) 



3. The projection of il 

As was demonstrated by Ritzwoller & Lavely (1991) and 
SCDT, this choice of orthogonal fitting polynomials V s 
leads to a s such that they correspond one-to-one to a 
projection of onto polynomials W s (i.e. Eq. (3.21) of 
SCDT) : 



2na n i 
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so that the latitudinal functions (t = 2s + 1) are : 



so that, if the fl s (r) are reconstructed from a sequence of 
1-dimensional inversions using the a-coefficients, f2(r, 9) 
can be reconstructed from : 



Q(r,0) = Y^VL s {r)W s {9) 



(22) 



Since the a-coefficients are independent the associated 
projection functions W s must be orthogonal to the latitu- 
dinal response functions associated with the a-coefficients. 
Combination of Eqs. (2), (19) and (22) yields the following 
orthogonality conditions : 
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/ dcostf Y, 72+iG f i n (0)WV(0) - gi,8 aa > 
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} 
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so that the kernels K^ s (r) are given by : 
K%\(r) = F?\r)g ls + F? l (r)g 2s 



(24) 



If conditions (23) are not satisfied there is 'cross-talk' be- 
tween a-coefficients for different s. Since the set W s is itself 
complete it follows from Eq. (23) that the linear sums of 
the G\ and of the Gi must be identical apart from a con- 
stant factor : 



jlm _/y(0 2s +l 



m— 1 
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Combination of Eqs. (9) and (25) yields the equation : 

i*Gf 

du 2 



(l-„»)ig_ = AG<') 



(26) 



Integrating the terms in Eq. (25) over u and using Eq. 
(10) shows that : 



A = -l or £ 7 ™ +1 =0 



The regular solution of Eq. (26) is : 



du 



(27) 



(28) 



for A = — 7}t(t + 1). It is convenient to choose A\ ' such 
that : 



(I) = (2t+l)vP 
' 2t(t + 1) 



(29) 



Gl 2(2 S + 2)(2 S + 1) (1 U) d,i^ +l( 



G« 2s+1 («) = -i(2 s + l)(2 s + 2)Gf )2s+1 ( 3 °) 

= -i(4s + 3)4^1(1 -" 2 )^2° s+ i(«) 

With the definition of Gi and G 2 through expressions (30) 
it is simple to show that : 



duG { l )2s+ \u) 



duG 



(l)2s+l. 



{u)=v < i] +1 5 sQ (31) 



By integrating all terms over u in Eq. (25), and by using 
Eqs. (10), (25) and (31), it follows that for all s ^ the 
coefficients 7 satisfy J2m=i Y2I+1 = as demanded by 
(27). The 

v 2s+i are determined by the normalization of 
the fitting functions V which is demonstrated in Sect. 4. 
Equations (30) can also be written as : 

Gf )2s+1 (u) = -l -ili±g^k- (l - u*)ipl +1 (u) 
1 v ; 2 (2s + 2)(2s + l) K ' 2s+lv ' (32) 

Gf 2s+1 («) = \vf s+1 {As + 3) (1 - u*) hpl +l ( U ) 

From the orthogonality properties of associated Legendrc 
polynomials (cf. Gradshteyn & Ryzhik, 1994) it follows 
that a natural choice for the projection polynomials W s , 
independent of I, for 17 is : 



W s {u) = -(l-u*)-^Pl +1 (u) 



(33) 



From Eq. (33) it follows that then conditions (23) are sat- 
isfied with : 

(0 

91 s = ^2 S +l 



92s = -^l +1 (2 S + l){2 S + 2) 



4. Normalization of the kernels 



(34) 



The freedom of choosing the factors i4«+i introduced in 
Eq. (29) is directly related to the freedom of choosing the 
normalization of the fitting functions V s . Combination of 
expressions (25) and (32) yields : 

1 (45 + 3)41+1 , , N i 1 
2 (2 fl + 2)(2 fl Tl) (1 - U ) ' P ^ (35) 



72 S +1 G 1 = - 



m— 1 



Multiplying Eq. (35) by W s i and integrating over u on 
both sides yields : 



E-.m aim 
728+1^25' + ! 



' r 

^72™+! duG[ m (u)WAu) 

m=l K ( 36 ) 
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where the factors /S^+i are introduced as shorthand no- 
tation for the integral. Using Eq. (20) it can be seen that : 



E <iWti = E K } + iM)' 



(37) 



m— 1 



From which it follows, using the orthogonality of the V 
polynomials, that : 



1 



+1 _ m 28+1 +1 
and therefore : 



^2i+i("»)=m J duG l l n (u)W s (i 



(38) 
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(l + l/2)x 
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Combination of the normalization condition T^s+iCO = ' 
with Eq. (39) yields : 



I = - 



l{2l + \) [(2Z-1)!!] 



2v. 



(0 

2s+l 



21) 



1 



Equation (40) is derived by using the explicit expression 
for P{(u) (cf. Gradshteyn & Ryzhik, 1994) : 



Pl(u) = (-l) l (2l -!)!!(! -u 2 ) 1 ' 2 



(41) 



The solution of the integral in Eq. (40) is (cf. Gradshteyn 
& Ryzhik, 1994) : 



J du(i-u 2 y- i / 2 Pi +1 ( U ) = 

2tiT(Z + 



(42) 



T(l + S + mi- S )T(s + l)T(-^-s) 



which means that the coefficients ?4s+i must satisfy the 
following expression : 

(I) _ ( n .j (2Z + l)!(2s + 2)!(Z + a + l)! 
2s+1 [ ' I s\{s + l)\{l-s-l)\{2l + 2s + 2)\ (43) 

for s < I — 1 

For s = Eq. (43) reduces to = 1. By using that 

Pj^u) = -(1 - m 2 ) 1/2 it is easy to show that v[ l \m) = m 
which is identical to the result obtained by SCDT in their 
Appendix A. 



5. Constructing the fitting functions 

In Appendix A of SCDT a relatively straightforward tech- 
nique is given for constructing the orthogonal fitting poly- 
nomials V . It is also possible to obtain the fitting func- 
tions by a different iteration scheme, taking Eq. (39) as 
the starting point and making use of one of the functional 
relationships between Legendre polynomials of varying or- 
der and degree (cf. Gradshteyn & Ryzhik, 1994) : 

p "-t 1 - P "+V = ( 2l + -u 2 *T ( 44 ) 

When this relationship is applied to Pi s+1 in Eq. (39), this 
yields : 

2s+i{ m ) — t>2,-dm)+ ^ ^ (i + \m\y: 



u 2s+l 



(1 + 1/2) f duPliu) [Pr(u)Y 



(45) 



-l 



The integral can be evaluated by expanding the P2 S in 
terms of power series of (1 — u 2 ) and making repeated 
use of Eq. (44), now applied to the P™. This leads to a 

weighted sum of a set of integrals of [P ; ™^' fc '] 2 where k 
runs from — s to +s. The integrals are easily evaluated 
from the normalization of associated Legendre polynomi- 
als. The evaluation of the weighting factors in terms of 
k, I, to, and s is rather cumbersome however, and there is 
a more convenient route to the same result. Integrals of 
products of three associated Legendre polynomials also oc- 
cur regularly in quantum mechanics when adding angular 
momenta (cf. Merzbacher, 1970). It is straightforward to 
demonstrate that the integral in Eq. (45) can be written in 
terms of Clebsch-Gordan or Wigner coefficients, the con- 
struction schemes of which can be found in the literature 
(cf. Abramowitz & Stegun) . Equation (45) then becomes : 



s(2l + 2s + 1) 
'(2s + !)(/-«) 



P£ ) -iM + 



m(2Z + l) 



■2.+i (46) 
(U00\U2s0)(llm -m\U2sO) 



where the angular brackets are the usual notation for the 
Clebsch-Gordan coefficients. 



6. Implications for SOLA inversion methods 

All linear methods to perform inversions for rotation as a 
function of both radius and latitude reduce in algorithmic 
form to solving a large set of linear equations, ie. inverting 
a large matrix. The method of Subtractive Optimally Lo- 
calized Averages (SOLA) (Pijpers & Thompson, 1994) is a 
possible inversion method which has proven to produce re- 
liable results in other helioseismic problems. In the SOLA 
method the (symmetric) matrix to be inverted has as its 
elements integrals over r and u of all the cross products 
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of the individual kernels. For the 2-D inversion for solar 
rotation the number of elements in this matrix is of the 
order of 20000 x 20000 or more, which is prohibitive for 
the direct use of SOLA on this problem. However, a direct 
consequence of obtaining the explicit expressions (32) for 
the latitudinal kernels is that it can be demonstrated that 
this matrix has a sparse character. This can be used to 
make the direct use of SOLA quite feasible. 

The matrix elements for the SOLA method are : 



3 a 



dr / du K n i s (r, u)K n > Us ,(r, u) 



(47) 



For convenience of notation the single multiplet index j is 
introduced, which runs over all combinations of the indices 
n, I. The K n i s can be written as : 

K nls (r,u) = F^(r)G^ 2s+1 (u)+F^ l (r)G^ 2s+1 (u) (48) 

Substituting Eq. (48) into Eq. (47) for s and for s' leads 
to a sum of four integrals, and in each of the terms the 
integration over u is the following : 



j d U (4,s + 3)(l- U 2 ) 1 / 2 P 2 1 s+1 (4 S ' + 3)(l- U 2 ) 1 / 2 P 2 1 s , +1 (49) 



with the indices p, q being either 1 or 2, the explicit ex- 
pression for the diagonal blocks is : 



rjs j's _ \ ~* j-D Ajj' 
p,q 

with the factors ^ equal to 



Sll — u 2s+l u 2s+l 



s(2s-l) 



+ 



- „,(0 „,(0 



r D = 

S12 — 



S21 — S12 



S22 



,(0 .,('') 



(4s + l)(2s + 2)(2s + l) 
(s + 2)(2s + 3) 
(4s + 5)(2s + 2)(2s + 1) 
s(2s-l) (s + 2)(2s + 3) 

~ (8s + 2) (8s + 10) 

s(s + l)(2s + l)(2s- 1) 

(87+2) + 
(s + 2)(s + l)(2s + 3)(2s + l) 



(52) 



(53) 



(8s + 10) 



The lower off-diagonal blocks are the transpose of the up- 
per off-diagonal blocks : C js j ' s - 1 = &' s ~ x js . Therefore 
only expressions for the upper off-diagonal blocks are given 
here : 



Equation (49) is obtained by making use of the expres- Qi» j'a-i _ Cpg^iq 
sions (32) for the G functions. Making use of Eq. (44) this p>g 
integral can be rewritten as : 

with the factors equal to : 



(54) 



1 






j d« [Pi Pl +2 ] [Pi, - Pl, +2 ] = 


— ± 

S 3 -l a' 


-2 (2s + 2)!" 
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4s + 1 (2s-2)!_ 


5ss> 


2 (2s + 2)! 
4s + 1 (2s- 2)! 


2 (2s + 4)!" 
+ 4s + 5 (2s)! 


3s+ls' 


-2 (2s + 4)! " 




_4s + 5 (2s)! 





(50) 
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S12 



S21 



S22 



.(0 

; 2s4 

,(0 



.(*') 



W 2;5+l U 2s-l"g s 



hs - x 8s + 2 



(I) (O (s + l)(2s + l) 

"2s+l v 2s-l 



(55) 



(0 „(0 

2s+l u 2s-l 



8s + 2 
s(s + l)(2s + l)(2s- 1) 
8s + 2 



The right-hand side of Eq. (50) is derived by making use 
of the orthogonality properties and normalization of asso- 
ciated Legendre polynomials (cf. Abramowitz & Stegun, 
1972), where the 8 is the Kronecker symbol. From expres- 
sion (50) it follows that the integral of the cross-product 
of two mode kernels is not equal to only if s' is equal 
to s — 1, s, or s + 1. The matrix s ^ s is therefore block 
tridiagonal. 

Introducing the notation : 



AJJ 
■^pq 



dr F; l F q 



n'V 



(51) 



Each block has a number of elements equal to N mu i x N mu i 
with N mu i the number of multiplcts in the mode-set. 
The number of such blocks is equal to the number of a- 
coefficients available, which for currently available mode- 
sets is around 20. When compared with a direct 2-D SOLA 
inversion this procedure reduces the computer storage re- 
quirements by a factor which is roughly equal to the num- 
ber of available a-coefficients (~ 20), and the CPU-time 
required for inversion of the matrix is reduced by a factor 
which is roughly equal to the square of that (~ 400). 
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7. Conclusions 

In this paper it is shown that the process of fitting 
Ritzwoller-Lavely type a-coefficients to ridges in the (I, v)- 
diagram leads to rotational mode kernels for the a- 
coefficients which can be expressed in closed form, as func- 
tions of associated Legendre polynomials. These explicit 
expressions for the kernels facilitate inversions for solar ro- 
tation as a function of radius and latitude, since it clearly 
is not necessary to apply the iterative procedure outlined 
in Appendix A of SCDT to the splittings mode kernels in 
order to obtain the a-coefheient kernels. For convenience 
the results that are relevant for 2-dimensional inversions 
for solar rotation from a-coefficients are collected below. 
The equivalent expression of Eq. (2) is : 
l l 



dr / duK nls (r,u)tt(r,u) 



27ra„z2s+i 



where the 'mode kernels' are now given by 



K nls (r,u) ee (r)G? )2a+1 (u) + Ff(r)G 
with : 

ee p(r)r 2 



-,(2)2s+l 



(56) 



(«)(57) 



/Inl 



pnl E 

and : 



p(r)r 2 [L-^ nl (r)] /I nl 



(58) 



G 



(1)28 + 1 



1 (4s + 3)< J 



+i 



2 (2s + 2)(2s + 1) 



G 



(l)2s+\ 



1 



^j +1 (4 S + 3)(l-^)*B 



1 

2s+l 



2s+li u ) 
(«) 



(59) 



The factors w^s+i are given by Eq. (43). The latitudinal 
functions satisfy a normalization slightly different from 
that expressed in condition (10) : 
l l 

-l -l 



J AuG ( (' 



1 

/ 



duGf )2s+1 (u) 



(60) 



With these equations a 2-D inverse problem is defined 
which can be solved using methods identical to those used 
to invert for the solar rotation rate from individual fre- 
quency splittings. 

Furthermore it is shown that SOLA inversions for 2- 
D rotation can be speeded up considerably because the 
matrix to be inverted in the SOLA method is shown to be 
symmetric block tridiagonal. 
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